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ABSTRACT 

Observations with the Short Wavelength Spectrometer (SWS) onboard the 
Infrared Space Observatory (ISO) have led to the first detection of the methyl 
radical CH 3 in the interstellar medium. The u 2 Q— branch at 16.5 /iin and the 
-R(O) line at 16.0 /im have been unambiguously detected toward the Galactic 
center SgrA*. The analysis of the measured bands gives a column density of 
(8.0±2.4) xlO 14 cm~ 2 and an excitation temperature of (17 ± 2) K. Gaseous 
CO at a similarly low excitation temperature and C 2 H 2 are detected for the 
same line of sight. Using constraints on the H 2 column density obtained from 
C 18 and visual extinction, the inferred CH 3 abundance is (1.3 +2 ^) x 10~ 8 . 
The chemically related CH 4 molecule is not detected, but the pure rotational 
lines of CH are seen with the Long Wavelength Spectrometer (LWS). The 
absolute abundances and the CH 3 /CH 4 and CH 3 /CH ratios are inconsistent 
with published pure gas-phase models of dense clouds. The data require a mix 
of diffuse and translucent clouds with different densities and extinctions, and/or 
the development of translucent models in which gas-grain chemistry, freeze-out 
and reactions of H with polycyclic aromatic hydrocarbons and solid aliphatic 
material are included. 

Subject headings: Line: identification — ISM: abundances — ISM: lines and 
bands — ISM: molecules — Galaxy: center — Infrared: ISM: lines and bands 
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1. Introduction 

The methyl radical CH 3 is an important intermediate product in the basic ion-molecule 
gas-phase chemistry networks in the interstellar medium driven by cosmic-ray ionization 
(Herbst & Klemperer 1973). Together with CH and CH 2 , it is produced by a series of 
reactions starting with C + H3 — > CH + + H or the radiative association of C + + H 2 — > 
CH^ + hi/, followed by hydrogen abstraction reactions and dissociative recombination. 
Alternatively, it can be produced by photodissociation of methane (CH 4 ). Subsequent 
reactions of C + with CH 3 form one of the most important steps in the formation of more 
complex hydrocarbons. 

The ion which initiates the chemistry in the cold gas has been detected only recently 
toward the Galactic center by Geballe et al. (1999) at a surprisingly high abundance. This 
unique line of sight turns out to be an extremely valuable environment to study abundances 
in the cold low density interstellar medium, since even minor species like CH 3 may be 
detectable due to its long path. 

Although CH 3 is a simple species, it is difficult to obtain accurate laboratory 
measurements of its molecular parameters since, as a radical, it recombines very fast with 
other particles in a gas. Herzberg (1961) and Herzberg & Shoosmith (1956) were the first 
to determine that the molecule is planar, but definite proof came only from measurements 
of the out-of-plane bending mode v 2 at 16 \xm. by Yamada et al. (1981). 

Observations of the 16-16.5 fim wavelength range are strongly hampered from the 
ground due to the Earth's atmosphere. The first detections of CH 3 in space have become 
possible only using the Infrared Space Observatory (ISO) (Kessler et al. 1996). Bezard et 
al. (1998, 1999) have recently detected CH 3 in the atmospheres of Saturn and Neptune 
respectively, but no previous searches for the molecule in interstellar space have been 
reported. 
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2. CH 3 spectroscopy 

Since the CH 3 radical is planar and symmetric, it does not have electric dipole allowed 
rotational lines which could be detected in the (sub-)millimeter wavelength range. The 
planar nature also implies that the symmetric stretch v\ is infrared inactive and the 
asymmetric stretch v 3 at 3.16 /im relatively weak. Indeed, the transition dipole moment of 
the V3 band is found to be a factor of three weaker than the out-of-plane bending mode z/ 2 
(Triggs et al. 1992, Amano et al. 1992). 

To calculate the u 2 spectrum, the term energies were taken from Yamada et al. (1981). 
The nuclear spin of the H-atoms can couple either to a quartet or doublet state, with 
nuclear spin statistical weights of 4 and 2, respectively. Because CH 3 follows Fermi-Dirac 
statistics, the K — 3, 6, 9, . . . levels are quartet states, and the other K-values doublet states. 
The strongest Q-branch has N = K and is located at 16.5 /iin; the strongest other feature 
is the ®R(0) line at 16.0 /xm thanks to its favorable Honl-London factor. The band strength 
of (2.5 ± 0.8) x 10~ 17 cm -1 (molecule cm~ 2 ) -1 was taken from Wormhoudt & McCurdy 
(1989). The calculation of the spectrum was performed as described in Helmich (1996). 

The shape of the spectrum is very sensitive to the excitation temperature (see Figure 
7.14 of Helmich 1996). Besides the strong Q Q-branch and Q R(0) lines, many more features 
become visible at excitation temperatures above 25 K, most notable the satellite ^Q-branch 
at 16.53 /iin and the -P(2) line at 17.60 /iin. 

3. Observations and Data Reduction 

Observations were carried out in the SWS grating mode AOT06 (de Graauw et al. 1996) 
at a spectral resolving power R = A/AA « 1500 — 2200. The spectral range covering the 
CH 3 Q-branch and the -P(2) line has been measured on 1997 February 21 15:24:27-19:03:01 
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UT, whereas that covering the R(0) and R(l) lines has been obtained on the same day at 
19:03:45-20:55:41 UT. The SWS aperture size was 14" x 27" and has been centered on the 
position of Sgr A* RA 17 A 45 m 40 s .O, Dec -29°00'28".6 (J2000 coordinates), with the long 
side of the slit oriented within 1 degree of the north-south direction. Due to the rather large 
aperture size, the Galactic center sources IRS 1, 2, 3 and 7 also fall inside the beam (see, 
e.g., Geballe et al. 1989), whereas IRS 5 and 6 are positioned just outside the slit. 

Data were processed within the SWS interactive analysis system, based on the standard 
ISO pipeline OLP V8.7 products. The data reduction adhered to the recommendations of 
Salama et al. (1997). Raw data were rebinned to R = 5000, a value significantly larger 
than the actual spectral resolving power of the SWS to avoid losing spectral detail when 
convolving the observed data samples by the bin. 

The absolute calibration of the SWS data has about ±20% uncertainty on average 
longwards of ~15 /im (Salama et al. 1997). However, since our analysis is entirely based 
on spectra where the continuum is divided out, the actual uncertainty in the results is 
determined by the noise in the data rather than the actual calibration uncertainty. The 
main limitation of our analysis originates from the ±30% uncertainty in the CH 3 u 2 band 
strength (Wormhoudt & McCurdy 1989, Yamada & Hirota 1983). 

4. CH 3 results 

As shown in Fig. 1, the Q-branch at 16.5 fim and the R(0) line at 16.0 fim are clearly 
detected. This represents the first unambiguous detection of CH 3 in the interstellar medium. 
The upper limit on the -P(2) line at 17.60 /xm provides an important constraint on the 
temperature. Due to a blend with the [Nelll] 15.555 /j,m atomic fine structure line at the 
SWS grating resolution, no information from the -R(l) line at 15.54 fim could be obtained. 
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Both the Q-branch and the R(0) lines are shifted with respect to their expected LSR 
wavelengths by about —20 km s _1 . Although such a shift is close to the SWS wavelength 
calibration accuracy (Valentijn et al. 1996), a Vlsr = —30 km s -1 component of cold 
molecular gas has been reported previously by several authors from observations at radio 
and millimeter wavelengths with similar beam sizes (Serabyn et al. 1986 and Sutton et al. 
1990, CO; Pauls et al. 1996, H 2 CO; Serabyn & Giisten 1986, NH 3 ; Marr et al. 1992, HCO+; 
Giisten et al. 1987, HCN; Bolton et al. 1964, OH). In all cases several velocity components 
at ~ —50, —30 and km s _1 have been observed at much higher spectral resolutions. The 
relative strengths of these three components vary between the observed species with the 
km s" 1 component often the largest. At the SWS spectral resolution of ~ 150 km s~ x it is 
not possible to distinguish between these different velocity components, but our observed 
shift is consistent with a mix of them. The location of the absorbing gas can therefore not 
be attributed to one particular feature, but is possibly spread along the line of sight toward 
Sgr A* among spiral arms and molecular clouds. 

Fits of synthetic spectra to the data, as described in §2 and matching the SWS 
resolution, were performed for different excitation temperatures (10 to 50 K), different 
Doppler parameters (b between 1.5 and 30 km s _1 ) and column densities. The best fit to the 
individual absorption depths (Fig. 1 a,b) as well as their ratio (Fig. 1 c) is given in Table 1. 
The absorption depths are almost independent of Doppler parameter, and are mainly a 
function of column density. The ratio of the depths is a strong function of the excitation 
temperature. The inferred low excitation temperature of (17 ± 2) K from the 16.0/16.5 
fim ratio of (0.8±0.15) is consistent with the non-detection of the P(2) line. Because CH 3 
has no dipole moment, the populations of the lowest rotational levels are controlled by 
collisions, so that the excitation temperature is close to the kinetic temperature. 

The H 2 column density along the line of sight has been constrained by several sets 
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of observations. First, the measured extinction of 31 mag (Rieke et al. 1989) implies 
7V(H 2 )~ 2 x 10 22 cm -2 using N H /A V = 1.9 x 10 21 cm~ 2 mag -1 and assuming that at least 
half of the hydrogen is in molecular form. Second, the detection of at least one optically 
thin C 18 line (R(0) at 4.7716 /im) in our ISO-SWS observation together with the measured 
CO excitation temperature of 8-13 K implies iV(C 18 0) = (2 ± 0.5) x 10 16 cm- 2 . This is in 
good agreement with the analysis by Moneti & Cernicharo (2000) on the same data. Using 
16 0/ 18 0= 300 (Wilson & Rood 1994) and CO/H 2 =10- 4 implies iV(H 2 )^ 6 x 10 22 cm" 2 . 
We adopt A^(H 2 ) = (6 ± 3) x 10 22 cm -2 , leading to a CH 3 abundance with respect to H 2 of 
x(CH 3 ) = (1.3^ 2 ' 2 ) x 10~ 8 . Note that all abundances would be increased by a factor of 3 if 
the lower H 2 column density derived from the extinction is used. 



5. Related species: CH 4 , C 2 H 2 and CH 

The availability of the full SWS scan allows searches for other chemically related 
molecules. Specifically, the z/ 2 /z/ 4 dyad of CH 4 occurs around 7.7 fim and has been observed 
with the ISO-SWS toward massive protostars by Boogert et al. (1998). Toward the Galactic 
center, however, gas-phase CH 4 is not detected. Adopting the same excitation temperature 
as found for CH 3 , an upper limit for its abundance of <lxl0 15 cm -2 is found. Solid 
CH 4 is clearly detected by Chiar et al. (2000) toward Sgr A* with a column density of 
(3.0 ± 0.7) xlO 16 cm~ 2 . Thus, most of the CH 4 is in solid form, consistent with the low 
temperature. 

Detection of a blend of the pure rotational lines of CH at 149.09 and 149.39 /zm towards 
the Galactic center has been reported by White et al. (1999). We have re-analyzed the LWS 
observations carried out on 1998 February 20 10:11:34-11:06:44 in the LWS01 grating mode 
(Clegg et al. 1996) at a resolution of ~ 1500 km s _1 (Fig. 2 a). The LWS data reduction 
has been based on OLP V8.7 products and has been carried out within ISAP (Sturm et al. 
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1997). Outliers have been removed by iterative sigma clipping and the different scans have 
been flatfielded to their mean value by applying a second order polynomial offset to each 
individual scan. The fringing in the LWS data, present in all extended source observations 
has been removed by the dedicated module within ISAP. The inferred equivalent width for 
the unresolved doublet is (85 ± 5) km s^ 1 . Using the formulae from Stacey et al. (1987) and 
assuming T cx =17 K, this leads to iV(CH)^ (1.1 ± 0.1) x 10 15 cm" 2 . 

Finally, the z/ 5 band of gas-phase C 2 H 2 at 13.7 /im is clearly detected (Fig. 2 b). 
Following the analysis of Lahuis & van Dishoeck (2000), we find iV(C 2 H 2 ) = (5.5 ±0.8) x 10 14 
cm -2 assuming T ex = 17 K. Table 1 summarizes the results obtained from ISO observations. 
Note that the relative abundances of the molecules have much smaller error bars than the 
absolute values since the uncertainty in the H 2 column density cancels out. 

6. Chemistry 

The absolute and relative abundances of the observed molecules have been compared 
with a wide variety of models, including time- and depth-dependent models. None of the 
published pure gas-phase dense cloud models can reproduce the observations of all species 
(e.g., Millar et al. 1997, Lee et al. 1996, see Table 2). In general, the model CH 3 abundances 
are too low and the CH 4 abundances too high. Also, the model abundance of C 2 H 2 is 
significantly smaller than that of CH 4 , in contrast with observations. The only models which 
come close to matching the absolute and relative abundances of CH 3 and CH 3 /CH and 
CH3/CH4 are low-density translucent cloud models with n(H 2 )Ri 10 3 cm -3 and Ay s^few 
mag, so that photodissociation of CH 4 to CH 3 and CH 2 is still efficient. Table 2 lists recent 
model calculations by Terzieva & Herbst (priv. communication) and results based on the 
models by van Dishoeck & Black (1986) and Jansen et al. (1995) using updated branching 
ratios for the dissociative recombination of the hydrocarbon ions (Andersen et al. 2000). 
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Low metal abundances are favored, to prevent destruction of the hydrocarbons by oxygen 
and by sulfur atoms and ions. Note, however, that even though the abundances at Ay 3 
may match the data within a factor of a few, the CH 3 column density in such models is 
only 1 x 10 13 cm -2 , nearly two orders of magnitude below observations. At the same time, 
the CH column density of 1.4 x 10 14 cm -2 is a factor of 10 below observations. Because of 
the strong depth dependence of the CH, CH 3 and CH 4 abundances, it is not possible to 
reproduce the column density ratios with these same models. The large observed column 
density suggests that there are several clouds along the line of sight. Some combination of 
low-density diffuse clouds to produce the CH and denser clouds to account for the solid 
CH 4 may explain those data, but the mix would have to be tailored very specifically to 
simultaneously approach the large column densities of CH 3 and C2H2. 

An alternative suggestion is to invoke turbulent chemistry, in which a high CH + 
abundance leads to enhancements of other hydrocarbons by 1-2 orders of magnitude (e.g. 
Hogerheijde et al. 1995, Joulain et al. 1998). However, the relative ratios of CH and CH 3 
are unlikely to change in such models. 

Given the detection of solid CH 4 and the low inferred temperatures, it is plausible 
that gas-solid interactions and grain-surface chemistry also play a role in producing the 
hydrocarbons. In this respect, the situation for CH 3 may be similar to that for NH in 
diffuse clouds (Mann & Williams 1984, van Dishoeck 1998). Conversion of atomic carbon to 
small hydrides on grain surfaces may be significant, but no model results exist yet for these 
conditions. Such models should also explain the C 2 H 2 abundances and lack of complete 
C 2 H 2 freeze-out. Alternatively, reactions of atomic H with PAHs and solid aliphatic 
hydrocarbon material, known to be present toward Sgr A* from the 3.4 /im absorption 
feature, may lead to CH 3 . Shock chemistry is not likely to be important for this line of 
sight because of the low temperatures. 
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Future high spectral resolution observations of CH 3 toward Sgr A* to constrain the 
velocity structure, as well as observations of CH 3 and other molecules in different types of 
diffuse clouds are needed to constrain the basic hydrocarbon chemistry. 

The authors are grateful to the SWS instrument teams, to E. Herbst, R. Terzieva and 
D.J. Jansen for updated model results on CH 3 and to W. Duley for inspiring discussions. 
This work was supported by DARA grants no 50 QI9402 3 and 50 QI 8610 8 and by 
NWO grant 614.41.003. CMW acknowledges receipt of an ARC Australian Postdoctoral 
Fellowship. 
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Fig. 1. — (a),(b): Data and synthetic spectra of the CH 3 observations; (c): Determination 
of the CH 3 excitation temperature T cx from the integrated absorption ratio of the 16.0 /im 
-R(O) line and 16.5 /im Q-branch. 

Fig. 2. — ISO spectra of Sgr A*: (a) LWS spectrum of the CH pure rotational line doublet; 
(b) SWS spectrum around the C2H2 1/5 Q-branch and the synthetic spectrum. 
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Table 1: 


Summary of results 






Molecule Col. Dens. 


T 

x ex 


Obs. Mode a 




[cm 2 ] 


[K] 




CH 


(l.l±0.1)xl0 15 




LWS01 


CH 3 


(8.0±2.4)xl0 14 


(17±2) 


SWS06 


CH 4 


<lxl0 15 


17 b 


SWS06 


C2H2 


(5.5±0.8)xl0 14 


17 b 


SWS06 


c 18 o 


(2.0±0.5)xl0 16 


8-13 


SWS06 



a For details of the observing modes see de Graauw et al. (1996) and Clegg et al. (1996) 
b Assumed excitation temperature 
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Table 2: Comparison of observed abundances with models" 



Model 


CH 


CH 3 


CH 4 


C2H2 


Obs. 


18 + 2 7 2 (-9) 


13 + 7 (-9) 


<17(-9) 


9 ^(-9) 


LBH96 b 


0.26(-9) 


0.20(-9) 


180(-9) 


15(-9) 


TH00 c 


44(-9) 


8.6(-9) 


5.6(-9) 


9.2(-9) 


VDB00 d 


3.2(-9) 


6.4(-9) 


18(-9) 


2.0(-9) 



a All abundances with respect to H 2 

b Lee et al. (1996) new standard model with n H = 10 4 cm -3 and T=10 K at steady state 
with low metal abundances 

c Terzieva & Herbst (priv. comm.) translucent cloud model with n H = 2 x 10 3 cm~ 3 , T=10 
K and Ay=S mag at steady state with low metal abundances 

d Updated models of van Dishoeck & Black (1986) and Jansen et al. (1995) for % = 2 x 10 3 
cm -3 and A v =3 mag at steady state with low metal abundances 



